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Abstract 

We analyse pion-nucleus reactions in a microscopic transport model of the 
BUU type, which propagates nucleons, pions, deltas and 1440) -resonances 
explicitly in space and time. In particular we examine pion absorption and 
inelastic scattering cross sections for pion kinetic energies T^^ = 85 — 315MeV 
and various target masses. In general, the mass-dependence of the experimen- 
tal data is well described for energies up to the A-resonance (~ 160MeV) 
while the absorption cross sections are somewhat overestimated for the higher 
energies. In addition we study the possible dynamical effects of delta- and 
pion-potentials in the medium on various observables as well as alternative 
models for the in-medium A-width. 
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1 INTRODUCTION 

Pion production in nucleus-nucleus reactions has been proposed in the early 80's as an 
observable to test the nuclear equation of state H, ^ |[ . Since then, and also for more 
general questions about the properties of hadronic matter at high density, a lot of 
experimental and theoretical studies on pion production in heavy-ion collisions have 
been performed 0, |^, |, |, |[. Nowadays it is widely accepted that pion production 
is not the most suitable process for studying the nuclear equation of state (EOS), 
especially in inclusive experiments, since the pions, which are created in the hot 
and dense reaction zone, are absorbed and reemitted several times before they leave 
the system such that the final pions emerge from densities zones p ~ 0.5po 0- 
Nevertheless, pion production is still an important field of present nuclear research 
due to the pion's important role in the overall reaction dynamics and the equilibration 
of the system. 

The most elaborate theoretical approaches for the description of pion production 
are microscopical kinetic models which include the propagation of pions and nucleon- 
resonances as well as their mutual interactions p], y]. Basically there are three differ- 
ent microscopical realizations: the Intra Nuclear Cascade (INC) [|T0| , [TI] , [13| , Quan- 



tum Molecular Dynamics (QMD) ||14|, |T^, and the Boltzman-Uehling-Uhlenbeck 



model (BUU)lTg, 0, 0, ig In the INC model the nucleus-nucleus collision is 
simulated as the sum of all individual nucleon-nucleon collisions without taking into 
account self consistent mean-field potentials and Pauli blocking for the collisions. 
The QMD follows the same scheme as the INC, but takes into account the Pauli 
blocking in the collisions and a nucleus potential which is calculated as the sum of 
all two-body potentials. In this paper we will use the BUU-model which is quite 
successful in describing the experimental data on pion production in proton-nucleus 
as well as nucleus-nucleus collisions. Since the pion spectra and pion yields turn 
out to be quite insensitive to details of the treatment of the vrA^A-dynamics in the 
BUU codes for heavy-ion reactions, we study in this paper explicitly pion-nucleus 
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reactions. We note that similar studies have aheady been performed with the INC 



model several years ago pTl E2l [13 



From the pion physics point of view there has been a vivid interest in describ- 
ing pion-nucleus reactions, because pion-nucleus reactions are expected to provide 
information on the properties of the strong interaction and excited states of the nu- 
cleon in the medium. Two major aspects are important for the understanding of 
pion-nucleus reactions: first, nuclear structure effects are important when looking 
for exclusive processes, i.e. detecting also the final state of the nucleus, and, second, 
the reaction mechanism itself which determines also inclusive reactions. Especially 



the question of the absorption mechanism is widely discussed [0, |2^, ^ ^ in con- 
nection with multi-nucleon contributions and the properties of the delta resonance 
in the nuclear medium. 

The most successful model up to date to describe inclusive pion-nucleus data is 
the model from Salcedo et al. In a first step these authors calculate micro- 

scopically density-dependent probabilities for pion quasielastic scattering and pion 
absorption in the nucleus, also taking into account three-body processes. With these 
probabilities they perform a pion cascade calculation using nuclear density profiles 
without, however, propagating the nucleons (or A's) explicitly. We will compare our 
results with their calculations. 

In this paper we have the twofold goal to test the validity of the pion and res- 
onance dynamics of our BUU model and to investigate the reaction mechanism 
in pion-nucleus reactions. In this respect we compare our calculations with total 
and differential cross sections for inclusive measurements in the pion energy region 

= 85MeV — 315MeV for pion absorptive and inelastic processes. This energy 
regime is of particular interest since here the delta is predominantly excited in the 
nucleus. To address the above mentioned questions, we proceed in the following way: 
in section 2 we review the basic inputs of the BUU model, especially those which are 
relevant for pion-nucleus dynamics. In section 3 we discuss the results of our model 
in comparison with various experimental data, compare with previous calculations 
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and study the dependence of observables on the input of our model. In section 4 we 
give a summary. 



2 THE EXTENDED BUU MODEL 

2.1 Basic equations 

The dynamical description of hadron-nucleus or nucleus-nucleus reactions is based 
on the equation of motion for the time evolution of the nucleon one-body phase-space 
distribution /(r, 

df{f,p^,t) ^ _ ^^^^^.^^^^^^^^ .^.^^ ^ I[f{f,p,;t)] (1) 



dt m 



with 



4 /" r da 

I[fir,pi;t)] = J d^p2d^P3 J ^^4^12^ 

SHpi + P2 - Ps - P,) X (/3/4/1/2 - /1/2/3/4) (2) 

where U{f) is the nucleon mean-field potential, da/dQ is the differential nucleon- 
nucleon cross section, /j = 1 — f{f,pi]t) the Pauli blocking factors and vu the 
relative velocity between nucleon 1 and 2. Eq. (|1]) is called the Boltzmann- 
Uehling-Uhlenbeck equation (also known as Vlasov-Uehling-Uhlenbeck, Boltzmann- 
Nordheim or Landau- Vlasov equation). Its static homogeneus solution is the Fermi- 
Dirac distribution function. Details of the derivation of eq. (|I|) can be found in ref. 
p8| , Physically the Ihs of eq. (|l]) represents a Vlasov equation for a gas of nucle- 



ons moving independently in the mean field U. The rhs of eq. (0) is the two-body 
collision integral /[/] which describes changes of the phase-space distribution due to 
nucleon-nucleon collisions incorporating Pauli-blocking for the final nucleon states. 

Above the pion threshold inelastic processes become more and more important in 
nucleon-nucleon collisions such that pions, A's and even higher resonances have to 



be included in the model. Therefore, in 16, 17] we have extended the BUU model to 



4 



coupled transport equations for pions, etas, A-, A^(1440)- and A^(1535)-resonance 
distribution functions [I^ (see also [0). Since for the present study only nucle- 
ons, deltas, A^(1440)'s and pions play a role, we discard the explicit propagation of 
etas and higher resonances. 

Schematically one can write down the coupled equations used in our model in the 
following way: 

nf — -L -L -L J- 

DfA = I^a + Ia^n. (3) 

where Df always denotes the differential operator for the Vlasov equation and the 
collision term on the rhs of eq. (Q) incoporates the following processes: 

Ina^ In a NA ^ NA and AN ^ NN 

/^Tv-^iviv- NN* ^ NN* and A^A^* ^ A^A^ 

(4) 

^A^Ntt^ ^A'^Ntti ^A^N-k ^ ^ 

denoting A^(1440) by A^*. 

In the eqs. (0) all direct processes for pion production and pion absorption are 
neglected. Since it is known that about 20% of the pions in nucleon-nucleon collisions, 
especially at low energies, are produced in direct processes, we have studied the effect 
of the direct channel by taking into account also the collision terms 

^NN^NNtt l!^NN-^NN (5) 

for the change of the nucleon distribution function and the terms 

^NN-^NNn I^NN^NN (6) 

for the change of the pion distribution function due to collisions. 



Since in pion-nucleus reactions in the energy regime considered single resonance 
excitations clearly dominate, we neglect resonance-resonance scattering which, how- 
ever, is included for the nucleus-nucleus case |l3T |. 



2.2 The test— particle method 

The standard way to solve the coupled nonlinear integro-differential equations is to 
discretize the distribution function, i. e. to substitute the continuous distribution 
functions by a finite number of test-particles, i. e. 

1 N-A 

fir,p;t) = -J25ir-m)S{p-m) (7) 

i 

where denotes the number of test-particles per nucleon and A denotes the total 
number of nucleons in the reaction. The ansatz (5) is a solution of the Vlasov equation 
if fi{t) and Pi{t) follow the classical Hamilton-equations, 

dt ' 

§ - ^ • (8) 

dt rrii 

Thus the problem reduces to a classical time evolution of a system of a finite 
number of test-particles. By discretizing the time t we solve this simplified problem 
as follows: For given positions and momenta of the test-particles at time t their 
evolution aX t + 5t reads: 

p(t + 6t) = p{t) - 6t Vr=-f/(r, t) , (9) 
f{t + 6t) = m + 6t(^P^^ . (10) 

Instead of p/m we use p/E in the Hamilton-equation as appropriate for relativistic 
kinematics. We calculate the effect of collisions, resonance decays or absorptions 
between t + 6t/2 and t + 3 6t/2 by assuming straight line trajectories during this 
time interval. For the simulation of inelastic reaction mechanisms we use the parallel 
ensemble algorithm (ref. |2^) in which collisions are only allowed between particles 



in the same ensemble. We note that the isospin degrees of freedom are accounted for 
by separate distribution functions, but Coulomb-effects are not taken into account. 



2.3 Mean— field potentials 

The particles inside the nuclear medium propagate in different mean-field poten- 
tials as described in eqs. (P1),(|9|), ([T0|) . The nucleon mean- field potential is evaluated 
selfconsistently on the basis of a density-dependent Skyrme potential as: 

?7;v(p) =0.75toP + t3P" (11) 

with to = -1177.8Me\^//m^ = 1845.5Mey//m^ and r = 5/3, which leads to a 
binding energy of —IGMeV per nucleon at po ~ and an incompressibility 

K = 308AMeV. These parameters give quite realistic, stable density distributions 
for finite nuclei in the surface region. 

The delta and A^(1440) resonance are in general assumed to propagate within the 
same mean field as the nucleons if not explicitly noted otherwise. Deviations from this 
average potential will not only influence the propagation of the resonances R through 
the nucleus, but also change the collision probabilities; e.g. a difference between the 
nucleon and the resonance potentials will shift the energy in the RN — > A^A^ reaction 
and the peak position in the ttN —>■ R process. Thus to exploit the dynamical 
influence of an alternative delta potential we also incorporate different prescriptions 
for this mean field and discuss the implications in comparison to the experimental 
data. 

In our simulations we use two approximations for the delta potential. The first 
model relates to elastic pion-scattering experiments where one finds a delta potential 
Va ~ —SOMeV at po Thus we have employed a delta potential which is just 



the nucleon potential shifted to —30MeV at normal nuclear density, 

Ua{p) = Om (GeV) . (12) 

I Un{Po) I 
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As a second model we used a delta potential as extracted by Ehehalt et al. in 
33| from the delta-hole model. A suitable parametrization of this delta potential is 



given by |^ 



U^{p) = -0.7p+1.75p^/\GeV) . (13) 

This potential is also —30MeV at po, but shows a different density dependence. 

Furthermore, in line with the delta-hole model also the pion propagates in its 
own potential field in the medium. We thus consequently also use the respective 



pion potential as evaluated and parametrized in our former work , where we have 
investigated the efffect of pion and delta potentials on asymptotic pion spectra and 
dilepton yields in heavy-ion coUisons. 

2.4 Simulation of the collision integrals 

The test particles in our approach collide with each other as in conventional cascade 
simulations with collision rates as described by eq. (^). The reaction probabilities 
are calculated on the basis of free cross sections for the different processes and im- 
plemented as explained in detail in [|l^ taking into account the Pauli blocking of the 
nucleons in the final state with proper isospin factors. The time-step used in the 
calculations is At = 0.5/m/c; we did not find any significant change in the results 
when decreasing At further. 

2.5 ttN cross section 

For the ttN —>■ A cross section we employ a Breit-Wigner formula 

_^maJ<lry T{M) 
3 , , / / \ \ 2 



with the momentum-dependent width 

Or) M \v 



r(M)= W ^ 44) r., (15) 
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where M is the actual delta mass, Ma is the peak delta rest mass (1232Mel^), q is 
the pion momentum in the rest frame of the delta, Qr the value of q for a delta mass 
Ma and = llOMeV. The function 

with (3 = SOOMeV cuts the width at high momenta. 



From the particle data table [35| we get 



= 200m6 

<-7 = 70mb (17) 
and based on isospin invariance of the strong interaction we use 

aX = 70m6 (18) 

in eq. (|l^) for the respective cross sections. 

The production of the A^(1440) is only possible in the Pu channel. In order to 
reproduce the vr+n cross section around the A^(1440) resonance we use the following 
form of the cross section for the nN — > A^(1440) process 

with cr^^^. = 30m6 and Tn* = 200MeV which yields the cross section (dashed 
line) in fig. 1. The sum of A and A^(1440) cross sections (eq. (|1^ with eq. (p!8D and eq. 
(^)) then gives the correct value of the total cross section (cf. fig. 1 (dashed-dotted 
line)). 

If a A or A^(1440) is created in a vrA^ collision, its mass is determined by the 
invariant mass y/s of the ttN system; if it is created in a A^A^-collision, then its 
mass is randomly chosen weighted with the proper Breit-Wigner distribution (see 
eq. (|29D). In our approach the resonances are treated as on-shell particles, however, 
with dynamically determined mass and lifetime (see below). 



2.6 Decay of resonances 

The decay of a resonance R — > t^N in free space is determined by its life time 
r(M)~^. However, inside the nucleus the resonance decay may be forbidden by the 
Pauli-principle. In addition, a resonance can also decay by the process -RA^ NN 
due to the presence of other nucleons. We will discuss the latter channel in the next 
subsection and describe here the numerical implementation of the 'free' decay inside 
the nucleus. 

Both resonances, A and A^(1440), decay with 99.5% and 65%, respectively, into 
a vrA^ final state. Thus in every time-step we calculate the decay probability into a 
TriV state in the rest frame of the resonance by assuming an exponential decay law 

Pdecay 6 ^ (20) 

where T{M) is the width of the resonance and At is the time-step in our calculation. 
By a Monte Carlo method we then decide, if the resonance actually decays in a given 
time-step. In case of a A we use eq. (^) for the width, for the A^(1440) a constant 
width r = 200Mey. Before we let the resonance decay we check if the outgoing 
nucleon is Pauli blocked by counting the number of testparticles in the final phase 
space cell as described in p8| . 

The A decay, furthermore, is assumed to be isotropic in its rest frame and the 
pion and nucleon momenta are determined by energy and momentum conservation. 
Nevertheless, to test the influence of an anisotropic delta decay we also have per- 
formed calculations with a p-wave angular distribution for the differential vrA^ cross 
section 

^~(l + 3cos'^e) (21) 

where G is the angle between the initial and final pion in the delta rest frame (c.f. 
section IsTSl) . 
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2.7 Direct pion emission and absorption 

Whereas in the VerWest-Arndt parametrization for the A^A^ NA —>■ NNn pro- 
cess c.f. section [27^ ) all pions are assumed to be produced via a delta or 
a A^(1440)-resonance up to ~ 2.5GeV, the pion may also be created directly 
in a NN collision. We calculate the direct pion production cross section, which is 
the production via an off-shell nucleon, in lowest order perturbation theory using a 
One-Boson-Exchange model for the nucleon-nucleon interaction. The full theory, 
that takes into account intermediate deltas and nucleons, reproduces the experimen- 



tal data rather well [0 . To extract the direct contribution to the production process 
we dropped all diagramms containing a delta in the intermediate state. (For a more 
detailed description of the calculations we refer the reader to ref. [^.) 

For the different isospin channels we use the isospin decompositon of the A^A^ — »• 
A^A^vr process (according to ref. p6[]): 



pp — i> ppn^ ail 
pp pnn+ 0-11 + o"io 

np — > rap7r°+ |(crio + aoi) (22) 

np nmT'^+ |(crii + aoi) 

np ppTi- + |(o-io + o-Qi ) 

where /, /' in ajj' denotes the initial and final isospin of the participating nucleon 
pair. By calculating the direct contribution to the pp — > pnn^, pp —>■ ppir^ and 
np —>■ npTT^ processes we determine all isospin channels using decomposition eq. 
(P^). Our result for the direct pp npir^ process is shown in fig. 2 (solid line). It 
has the following functional dependence on ^/s (in GeV): 



a, 



pp^npn' 



-(v^) = 26.2 - 40.6712v^+ 13.55566s . (23) 



For the other isospin channels we used the same energy dependence with the following 
ratios (calculated at y/s = 2.155GeV): 

'^pp—>npn+ ■ ^np^np-K^ ■ '^pp-^ppn^ 6.98 . 2.269 . 1 . (24) 
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Taking into account the direct pion production process, we simultaneously have 
to consider the inverse process nNN NN. We avoid the problem of calculating a 
cross section with a three-particle initial state by using the quasi 'deuteron' assump- 
tion; the pions are assumed to be absorbed on nucleon pairs as in the conventional 
picture for pion absorption ||25[. The 7r(2A^) NN reaction we then can calculate 
from the direct production cross sections by using detailed balance: 

(^7tNN-*nn - ■:^—'^nn~,nntt l^^J 

O Pi 



with 



pt = (s- (m,-2m„)')(s- (m^ + 2m„)')/(4s) (26) 
Pj = ' (27) 

which denote the initial (2A^) pair and final nucleon momenta in the center of mass 
frame. In fig. 2 we show the iT~^{pn) — * pp cross section used in the calculation 
by the dashed line. The isospin dependence discussed above leads to the following 
absorption ratios: 

-- 0.083 
-- 0.44 

-- 0.14 (28) 
-- 0.14 
-- 0.083 
-- 1 . 



^TT^ (nn)- 


■*np 


'^7r+ (np)- 


■*PP 




■*np 


(^TT+ (np)- 


-^pp 






^7r+(np)- 


-^pp 


^TT^(pp) — 


■*pp 


'^7r+(np)- 


-^pp 


^Tr~ (pp)- 


■^np 


C"7r+ (np)- 


■*pn 


^■K~ {np)— 


*nn 


^■K+{np)- 





2.8 NA ^ NN and iViV(1440) ^ NN reactions 

For the A^A^ NA and A^A^ A^A^(1440) processes we use the VerWest-Arndt 



parametrization of the cross section as described in detail in ref. |16]. The mass of 
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the resonance R, which is populated in our simulation in a A^A^-coUison, is chosen 
according to the probability function 

FiM^) = 1 ^^^(^) (29) 

^ ^ TT (M2 - Mi )2 + Mir(M)2 ^ 

with r(M) for deltas given in eq. (0) and r(M) = 200MeV for A^(1440). In a pion- 
nucleus reaction these resonance production processes are negligible, but the inverse 
processes A^A NN and A^(1440)A^ NN are very important in the context 
of pion absorption. In the literature the latter reactions are determined from the 
measured cross sections for the inverse channels by detailed balance using slightly 



different assumptions |0, |T^, that have to be discussed in more detail. The 

ambiguity arises because the measured cross sections given in ref. are averaged 
over the A-mass distribution. 

Application of detailed balance to a hypothetical delta with fixed mass leads to 

2| 



1 

for the A'^'^n — > pp cross section. In eq. (^) Pf,Pi are the final and initial momenta 
of the particles in the center of mass frame whereas the factor 1/4 is due to spin 
averaging and a symmetry factor for identical particles in the final state. 

Taking into account, however, that the delta has no fixed mass in our simulation. 



we get [17 



^nA++^pp - -J^f(m^)J_, d(cose)dM2^^'°'"^ • ^^^^ 

Since the VerWest-Arndt parametrization only gives the A-mass averaged total pion 
production cross section it cannot be used to calculate the inverse cross section in 
eq. (|3TD. Following our former work we therefore use 

4 pl j(v^-^)^ /(M2)rfM2 

with the A^A^ — * A^A cross sections from VerWest-Arndt. 
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To clarify the situation concerning the different expressions in the hterature for 
the A^A NN reaction we inserted in eq. ( pT] ) the result for the mass-dependent 
cross section dapp^nA++ / {dQdM'^) from a calculation of Schafer et al. I^^. The latter 
differential cross section for the pp nA"*""^ process is calculated on the basis of a 
OBE model using pion- and p-exchange which reproduces the experimental data for 
da/dM (cf. fig. 3) remarkable well. 



To investigate the validity of our model assumption in |jT^ we compare the results 
of eqs. (^),(^) and eq. (^Tj) calculated with the OBE model for the A+^n pp 
cross section for different delta masses M = 1180, 1232 and 1280Mel^ in fig. 4. Our 
cross section (solid line) eq. (0) is in good agreement with the OBE result (dotted 
line) for larger invariant energy ^/s, but overestimates the absorption cross section 
close to threshold. In fig. 4 also the cross section as proposed by Danielewicz and 
Bertsch in ref. |3^ is shown by the dashed-dotted line. It fits better at small ^/s but 
overestimates the absorption significantly at higher energy ^/s. The naive detailed 
balance ( pOf ) (dashed line) underestimates the cross section at almost all energies. 
For the dynamical implications of the different models in pion-nucleus reactions we 
refer the reader to sections ^?T] and 

Finally, using isospin factors and taking into account symmetry factors for iden- 
tical particles in the final state, we obtain the total cross sections of the other isospin 
channels as 



tot 
'^nA+- 


■^pn 


^ pA'^—^pn 


„tot 
2^nA++- 






tot 
^pA+- 


^pp ~ 


^tot _ 


^ ^iot 
2^nA++- 


-^pp 




tot 
^pA~- 




- 

'^nA++^pp 






(33) 



In case of the A^A^(1440) NN process we use eq. (^) with a constant width 
Fat. = 200MeV and the A^A^ A^A^(1440) cross section from the VerWest-Arndt 
parametrization. Additionally, we also check for Pauli blocking of the outgoing nu- 
cleons; however, due to the large momenta transfered to both nucleons the blocking 
effect is found to be of minor importance in pion-nucleus reactions. 
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2.9 NN, NA and iV(1440)iV elastic scattering 



For these processes we use the conventional Cugnon parameterization ITT], ^ 

35 

"^'^'^ = 1 + lOOVi'(GeV) '^^) 
with y/s' = v/i- 1.8993 for a NN colhsion and y/s' = y/s -0.938 - Mr for nucleon- 
resonance elastic scattering. The anisotropic angular distribution used is adopted 
from [jl6 . 



2.10 Inclusive cross sections 

To calculate inclusive cross sections for a pion-nucleus reaction we perform an explicit 
impact parameter integration, 

a = J 2Mh Nreac{h) , (35) 

where Nreacip) denote the impact-parameter dependent particle multiplicities as ob- 
tained from simulations with fixed h for the pions. 

This technique is identical to the one employed for proton-nucleus or nucleus- 
nucleus reactions so that we describe all reactions on the same footing. In ref. |1^, |3^ 
we have already shown that the experimental spectra on pion production in heavy- 
ion collisions can be reproduced very well, except for some small differences at small 
pion transverse momenta. If this also holds for pion-nucleus collisions will be studied 
in the next section. 



3 RESULTS FOR PION-NUCLEUS REAC- 
TIONS 

3.1 Total absorption cross sections 

As noted in the introduction the mechanism of pion absorption in nuclei is not 
well understood. Energy and momentum conservation rule out absorption of a pion 

15 



on a free single nucleon. However, due to the large momentum mismatch in the 
reaction the single-nucleon absorption is very unprobable in pion-nucleus reactions, 
too. The most important mechanism is thus the two-body absorption, where the 
pion is absorbed by a pair of nucleons. The quantitative contribution of higher order 
processes and especially that of three-body processes is still a matter of debate. 
For heavier nuclei the three-body process is controversely discussed to contribute 
between 10% and 50% of the total absorption cross section p6| . 



In the transport simulations the pion absorption mainly proceeds via a two-step 
mechanism: a pion first is captured in a A state and this A can be absorbed in a 
second step via the channel AA^ — A^A^. This two-step process clearly simulates 
a two-body absorption mechanism. However, by allowing the A to decay to t^N 
and the outgoing pion to be captured again in another A (subsequent two-body 
reaction chain), we incorporate also higher order absorption mechanisms with on- 
shell intermediate pions. 

In fig. 5 we show the calculated mass dependence of the pion absorption cross 
section at the pion energy T^- = IQbMeV in comparison to the experimental data. To 
illustrate the influence of the AA^ A^A^ reaction we have performed calculations 
with the different detailed balance prescriptions discussed in section |2.8| . We find 
that our previously used formula eq. (0) [|I^ (solid line), the OBE results eq. ( |3TD 



(dotted line) and the prescription from ref. |]5S| (dashed-dotted line) lead to the same 
accuracy in reproducing the experimental data, whereas the naive detailed balance 
(eq. (pOD) misses about 50% of the cross section. This reflects the fact, that the 
detailed balance prescription of eq. ( |5D| ) underestimates the elementary absorption 
cross section over the whole energy regime (see fig. 4) and explains why the cascade 
calculation of [0 underestimated the absorption cross section by a factor of 2. 

Turning now to the energy dependence of the absorption cross section (fig. 6), we 
find that we can describe the cross section at low energies up to the resonance region 
very well for all target masses. For higher energies (T^ = 315MeV"), however, our 
calculation (solid line) overestimates the data. This disagreement becomes larger for 
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heavier nuclei. It is worthwhile to note that even in this T^^ regime the A-dependence 
of the absorption cross section is well described. The energy dependence of the pion 
absorption is not altered when using the OBE model (dotted line) or the prescription 
from ref. ||3^ (dashed-dotted line) (see section |2.8| ) for the AA^ NN reaction. 

Compared to other calculations for the total absorption cross section as a function 
of target mass and kinetic energy of the pion ||2^, ^ our results are of comparable 
quality. The INC calculations have problems in reproducing the mass dependence of 



the total absorption cross section 123] and also the absolute value in the resonance re- 



gion, whereas the calculations of Salcedo et al. ||2^ also overestimate pion absorption 
at the higher pion energies for heavier nuclei. 

To simulate the effect of a direct three-body absorption process in our approach, 
we used a density dependent modification for the AiV NN cross section of the 
form: 

^'aN^NN = (1 + 3 — )crAAr^ArAr , (36) 

Po 

thus incorporating the density dependence of three-body pion absorption in an ad- 
hoc way. The result for the energy-dependent absorption cross section on and 
'^^^Bi with the density-dependent cross section ( PBD is shown in fig. 7 (dashed line). 
We find only a minor three-body effect for the "^^^Bi case. This is due to the large 
diameter of the i?i-nucleus; most deltas inside the i?z-nucleus scatter several times 
with nucleons and are absorbed. An increase of the AA^ A^A^ cross section as 



introduced in eq. (36) just leads to an earlier absorption and not to an increase of 
the total pion absorption cross section. 

For the ^^C nucleus we find a bigger effect of the three-body absorption process 
as simulated by prescription (^); fig. 7 (dashed line). In this lighter nucleus deltas 
do not scatter as often as in the Si-case and can leave the nucleus without being 
absorbed. Since the escape probability of the A's increases with their velocity, modi- 
fications of eq. (|3^) for the AA^ NN process affect the total pion absorption cross 
section dominantly at higher energies, thus overestimating the total pion absorption 
cross section even more. 
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Salcedo et al. [^] report a larger contribution of the three-body absorption pro- 
cesses of about 30%, but these calculations would lead to similar results as ours, if 
they omitted their three-body absorption probabilities. Similar to our case without 
three-body absorption probabilities the pions in their calculation would be absorbed 
further inside the nucleus via a two-body process. We thus conclude that the total 
inclusive pion absorption process in heavy nuclei is not sensitive to true many-body 
effects. It seems that subsequent two body steps, taking into account off-shell deltas, 
can simulate the many-body effects. 

We additionally note, that we have also studied the influence of the A^(1440) for 
the absorption process and have found only a 10% difference for T^^ ~ 300MeV when 
performing the calculations with and without the A^(1440) resonance. For smaller 
energies the contribution of the A^(1440) is negligible. 

As mentionend in section |2.7| we studied also the effect of direct (NN) pion pro- 
duction and absorption. We find a 10% enhancement for the absorption cross section 
at all pion energies in the ^^C case, whereas in the "^^^Bi case the pions are absorbed 
earlier in the reaction due to the incorporation of the direct process, but the total 
absorption cross section is not changed (see discussion on many-body absorption 
above). 

Summarizing this section we may state that we reproduce the experimental data 
on pion absorption quite well, for light nuclei at all energies and for heavy nuclei 
for energies up to the A-resonance region; for heavy nuclei above the A-resonance 
we obtain a too large cross section. The origin of this deficiency may be in the 
cascade description itself. We find that for central (tt. A) collisions the absorption 
cross section even has a minimum at T^r = IGOMeV because at this resonance energy 
the A-excitation cross section is very large and the delta can thus be formed already 
in the nuclear surface, where the probability for the decay pion to escape the nucleus 
again (inelastic scattering) is very large. The calculated energy dependence of the 
absorption cross section, which shows a different behaviour, is thus largely determined 
by the impact parameter dependence. 
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3.2 Total inelastic cross section 

Inelastic pion-nucleus scattering occurs in competition with the absorption process. 
Thus to understand the pion-nucleus reaction both processes must be described 
within the same approach. 

As shown in fig. 8 our calculations reproduce the mass dependence of the total 
inelastic cross section for pion energies [T^^ = lQ5MeV) (solid line). This also holds 
when employing the OBE model for the AA^ NN channel (dotted line) or the 
prescription by Danielewicz |^ (dashed-dotted line). 

When looking at the energy dependence of the inelastic cross sections, which is 
shown in fig. 9 for ^^C and "^^^Bi, we see that for the "^^^Bi case we underestimate 
the cross section by about a factor 2 at all pion energies; for this only happens 
for higher energies above resonance. 

A detailed phase-space analysis offers the following kinematic reason for these 
deficiencies: Generally most of the pions that are scattered inelastically emerge from 
the first A generation which is dominantly created in the surface of the nucleus; 
approximately 50% of the pions from the first A generation leave the system again 
and are not absorbed. When increasing the pion kinetic energy above the resonance 
region the cross section for delta formation becomes smaller and consequently the 
first delta generation is created further inside the nucleus, such that pions from its 
decay cannot leave the nucleus as easily as in the resonance region. Furthermore, 
these secondary pions are shifted down in energy (see section |33| ), which leads to an 
effective increase of the reabsorption cross section. As a consequence the secondary 
pions in the interior of the nucleus are almost completely absorbed, which leads to a 
reduced inelastic pion cross section and a correspondingly enhanced absorption cross 
section for higher incident pion energies. 

At low energies, on the other hand, the first generation of deltas has a larger life 



time due to their low mass (see section |3.5| ). This does not change the picture just 



described in heavy nuclei, but in light nuclei some of these deltas can thus decay 
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outside the nucleus and the secondary pions can escape again. 



3.3 Differential inelastic cross section 

Since in the resonance region the inelastic cross sections are well reproduced, we 
proceed with calculating also differential cross sections to get further information on 
the initial delta distributions. In fig. 10 the results of our simulations for 7r++^^C and 
^+_l_209 2j^ at = IGOMeV are shown for the differential inelastic cross section (solid 
line). Except for the forward angles the data for "^^^Bi are rather well reproduced both 
with an isotropic and an anisotropic delta decay (see discussion in section p.6|) . For 
the ^^C target, however, the anisotropic delta decay gives clearly a better description 
of the data at backward angles (dashed line in fig. 10). The deficiency in forward 
direction is due to the fact that the forward scattered pions will be captured in a 
delta again so that it becomes rather improbable to leave the nucleus in forward 
direction. Furthermore, the forward amplitude is also expected to be sensitive to 



coherent quantum scattering this effect is completely missing in our incoherent 
transport approach. 

Using the more realistic anisotropic delta decay instead of the isotropic one we find 
that this effects the total cross section for absorption and inelastic scattering by at 
most 5% such that the basic conclusions of the previous sections remain unchanged. 
Again our results are comparable with those of Salcedo et al. (ref. [^) and are 
quantitatively better than the results from INC calculations 



3.4 {7T~^,Np) reactions 

A further test for the absorption mechanism in our model is the comparison to 
more exclusive experiments where also the continuum protons are detected in a pion 



absorption event. In a recent study Ransome et al. present total cross sections for 
pion absorption on nuclei followed by 2p,3p emission for different energies and masses. 
We have calculated the (7r+, 2p) and (vr^, 3p) cross sections for = 150MeV, where 



20 



element 


(7r+, 2p&3p) 


(7r+, 2p&3j9) 


(7r+,3p) 


(7r+,3p) 




BUU 


exp. 


BUU 


exp. 




114 


110 


15 


14 


58 


272 


300 


30 


33 




333 


320 


46 


31 


208 


399 


400 


48 
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Table 1: Cross sections for (tt"*", Np) reactions in mh for T^^ = 150Mel^ from ref. [39]. 
We compare with the sum of 2p and 3p events, because these are the experimental 
data given under the heading total 2p. 

our calculation reproduces the total absorption cross section. For the continuum 
protons we applied the same cut in energy Ep > 23MeV as in the experiment. To 
calculate the cross section we proceeded as described in section |2.10| . For different 
impact parameters we counted the absorption events followed by two, three proton 
emission, respectively. 

Our calculated results are compared in table |l| with the data; they are in good 
agreement, especially for the lighter nuclei. For the heavier nuclei we get too many 
(7r+,3p) events, but for this channel the errors in the experimental data are bigger. 
In accordance with the experimental data we find that the contribution of (Tr"*", Np) 
processes with > 3 to the absorption is negligible. 

Thus we conclude that the average event characteristics are well described within 
our transport approach. 



3.5 Pion and delta dynamics 



In the calculations described so far the deltas have been propagated in the same 



potential as the nucleons. Recently [jT3|, |3^, there has been some discussion of 
the properties of the delta in the nuclear medium. In order to explore the physical 
consequences of any in-medium changes we have performed calculations changing 
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artificially the peak mass and width of the resonance. To be more precise, in the 
width r/^{M) we have changed the parameter in eq. (|15|), which consequently 
influences all other cross sections that depend on r(M). In varying systematically 
the peak mass of the delta and its width at maximum we find that both changes do 
not improve the energy-dependence of the absorption cross section discussed earlier 



in sections |3T^ and The best reproduction of the absorption data is achieved 
for the free width and the free delta peak mass (solid line in fig. 11). Thus our 
model analysis of pion absorption cross sections indicates no need for any in-medium 
changes of the A-resonance. 



In a recent publication Sneppen and Gaarde |^ have argued that the analysis 
of charge exchange reactions in terms of a cascade model requires an in-medium 
delta mass of Ma = 1202MeV' and a width of T{Ma) = 200MeV. In a calculation 
with these parameters we find for the ^^C target (fig. 12, (dashed line)) that the 
agreement with the experimental data at higher energies is better compared to the 
original calculation (solid line), but becomes worse for low T^. In addition, the 
inelastic scattering cross section at higher energies is not well described. For the 
'^^^Bi case (not explicitly displayed) there is not much change. 

Another in-medium effect for the deltas is the delta potential as already pointed 
out in section |2.3| . In order to explore the sensitivity of the pion-nucleus data to pion 
and delta selfenergies we have studied the influence of pion and delta potentials given 
in section ^]3| in our simulations. Whereas the effect of a pion potential changes the 
result only on a 5% level (not explicitly displayed), the pion absorption shows a larger 
sensitivity to a A-potential as shown in fig. 13. We find that the effective potential 
from the delta-hole model eq. ^ (dotted line) and the simple shifted potential eq. 
(0) (dashed line) lead to similar results. The calculated changes in the absorption 
cross section can be directly traced back to the energy shift in the nN reaction. 

Furthermore, to separate kinematical effects in the delta and pion dynamics in 
the pion-nucleus reactions we have calculated the differential A mass distribution 
for different A generations. The results in fig. 14 suggest the following picture: At 
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higher energies also high mass deltas are excited, as naively expected. These deltas 
in turn have a short life time (large T{M)) and decay fast. The pions from their 
decay are reabsorbed again and regenerate the A's although now at a smaller mass 
due to recoil effects. Thereby the A's cascade down in mass and loose about 70MeV 
in mass per generation. 

4 SUMMARY 

We have analysed pion-nucleus reactions within the extended BUU model that had 
originally been developed for the description of nucleus-nucleus collisions. Due to 
the pionic entrance channel we expected to be especially sensitive to the pion and 
delta dynamics employed and thus to be able to put a stringent constraint on these 
ingredients for the description of heavy-ion collisions. 

Apart from heavy-ion collisions we can reproduce the inclusive data for pion- 
nucleus reactions for a wide range of pion kinetic energies and target masses. Within 
the transport approach we have simulated pion absorption by sequential binary re- 
actions between pions, nucleons and delta's. The measured pion absorption cross 
sections are very well reproduced for pion energies up to the resonance energy; for 
higher energies we overestimate the absorption. The mass dependence, on the other 
side, is very well reproduced. The dominant part of the many-body absorption ef- 
fects seems to be simulated by the dynamics incorporated in our transport approach. 
The deficiencies of the model could essentially be traced back to the first generation 
of deltas. 

Apart from pion absorption we also reproduce the inelastic cross sections for 
light nuclei up to the resonance region quite well. For heavy nuclei we seem to 
underestimate the inelastic cross section at all energies. 

We obtain a good agreement for differential inelastic scattering and the exclusive 
(tt"*", Np) reactions in the delta resonance region which demonstrates that the overall 
event pattern is described quite well in this region. 
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In summary we find that the BUU transport approach developed for nucleus- 
nucleus collisions provides a rather reliable description of pion-nucleus reactions, 
too. The results obtained for the different channels in the pion-nucleus reactions 
are in their quality comparable to those of other recent model calculations for these 
reactions. 
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Figure Captions 



Fig. 1: Decomposition of the tt p cross section: Breit-Wigner cross section for 
A° (eq. dH)) (solid line) and Breit-Wi gner cross section for tt p A^(1440) 
(eq. (0)) (dashed line). Sum of eq. and eq. (|TP|) (dashed-dotted line). Experi- 



mental data are from ref. 41 



Fig. 2: Cross sections for the direct pp pmr^ process based on the OBE model 
(see text) (solid line) and the direct 7r+(]?n) — > pp process ( pSj ) (dashed line) as a 
function of -y/i. 

Fig. 3: The mass dependent cross section for pp A++n at y/s = 1.48 in the OBE 
model from S chafer et al. Experimental dat from ref. 



Fig. 4: The cross section for A'^^n pp within different descriptions for the 
detailed balance as a function on ^/s — M: naive detailed balance eq. (|30D (dashed 
line), calculation with the OBE model eq. ( PT| ) (dotted line), eq. ( P^ ) (sohd line), 
Danielewicz ||3^ (dashed-dotted line). The calculations are for delta masses M = 
1180, 1232 and 1280Me1/. 

Fig. 5: Mass dependent total absorption cross section for 7r+ at Tj^ = 165MeV 
within different descriptions for the detailed balance: naive detailed balance of eq. 
(0) (dashed line), calculation with the OBE model eq. (|3T|) (dotted line), eq. (|32D 
(solid line), Danielewicz ||3^ (dashed-dotted line). Experimental data are from ref. 



Fig. 6: Total absorption cross section for vr"*" on ^^C, ^^Fe and "^^^Bi as a function 
of the pion energy: Calculation with the OBE model eq. ( PTf ) (dotted line), eq. ( |5^ 
(solid line), Danielewicz ||3^ (dashed-dotted line). Experimental data are from ref. 
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Fig. 7: Total absorption cross section for 7r+ on ^^C and as function of the 

pion energy: Calculation with eq. (|3^ ) (solid line) and eq. (^) for the AA^ — > NN 
process (dashed line). Experimental data are from ref. 

Fig. 8: Mass dependent total inelastic cross section for 7r+ at T^^ = 165MeV: Cal- 
culation with the OBE model eq. (^) (dotted line), eq. (|3^ ) (solid line), Danielewicz 
Il38| (dashed-dotted line). Experimental data are from 



Fig. 9: Energy dependence of the inelastic absoroption cross section for 7r+ on 
and Calculation with the OBE model eq. (§T]) (dotted line), eq. (solid 

line), Danielewicz ||3^ (dashed-dotted line). Experimental data are from ref. ||43|| . 



Fig. 10: Differential inelastic cross section for 7r+ on (upper figure) and on 

^^C (lower figure) at = IGOMeV. Calculation with isotropic delta decay (solid 
line) and with anisotropic decay eq. (0) (dashed line). Experimental data are from 
ref. 



Fig. 11: Total absorption cross section for 7r+ on and ^^^Bi as a function of the 
pion energy. Calculations with different delta rest masses are shown in the left figure 
and with different Tj. in eq. (O) for all cross sections in the right figure. Experimental 



data are from ref. [43 



Fig. 12: Total absorption (upper figure) and total inelastic (lower figure) cross sec- 
tion for 7r+ on ^^C as a function of the pion energy. Calculations with the parameters 
from ref. [Q: delta rest mass Ma = 1202MeV and = 200MeV only in T{M) of 
eq. (pOl). Experimental data are from ref. ^3 . 



Fig. 13: Energy dependence of the total absorption cross section for 7r+ on 
and ^^^Bi. Calculations with different delta potentials: Va = Vn eq. (|lT]) (solid 
line), shifted delta potential eq. (0) (dashed line) and potential from ref. |Q eq. 
(O) (dotted line). Experimental data are from ref. [H3|. 
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Fig. 14: Delta mass distribution for different delta generations for tt^ on '^^^Bi 
at 85MeV (lower figure) and at 245MeV (upper figure). Free mass distribution eq. 
(^) (dashed line), first delta generation (solid line), second delta generation (daslied- 
dotted line), third delta generation (dotted line) and fourth delta generation (thick 
dashed line). 
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